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We review algorithmic methods for two-loop calculations in HQET, and the analogous 
methods for on-shell QCD, needed for matching HQET to QCD. 

1. Introduction 

In this talk, we review methods for two-loop calculations in heavy quark effective 
theoryS. In HQET, the heavy-quark propagator is S(k) = l/(fc • v + iO), and the 
quark-gluon vertex is igt a v^, with a heavy-quark velocity v — (1,0), in the rest 
frame. We shall discuss algorithmic methods, suitable for computer-algebra imple- 
mentation. All of them are based on dimensional regularization, with space-time 
dimension d = 4 — 2e, and integration by partsi, a powerful calculational method 
based on the simple observation that integrals of full derivatives vanish. Tensor in- 
tegrals are handled by decomposition into a suitable tensor basis, with coefficients 
that can be expressed in terms of scalar integrals, by solving linear systems of equa- 
tions. In all cases which we shall discuss, numerators can be expressed in terms of 
factors of the denominator. 

In Section 2, we recall the well-known integration-by-parts algorithm for calcu- 
lation of massless two-point integrals]. This example is used to introduce notations 
and to set up the pattern which will be repeated in other cases. In Section 3, we dis- 
cuss the corresponding method in HQETO. This method was generalized! to HQET 
three-point integrals, expanded in the velocity change v' — v, to any finite order, 
as discussed in Section 4. In order to obtain QCD/HQET matching conditionsEi, it 
is necessary to equate QCD and HQET on-shell matrix elements. The latter are 
trivial (apart from massive imark loops). The former require methods for on-shell 
massive QCD calculationsa-QElla, discussed in Section 5. 

2. Massless two point integrals 

We write the massless one-loop integral of Fig. [j]a, with arbitrary indices, as 

d d k 1 i(_p2)<*/2-e»i-a a 



{2n) d (-fc 2 ) Q i(-(fc+p) 2 )^ (47r) d / 2 



-G( ai ,a 2 ). (1) 
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It can be calculated by using the Feynman parameterization 
1 T(a + /3) [^(l-^da; 



a a bP T(a)T{(3) J [ax + 6(1 - x)] a +P 



(2) 



or by using Fourier transformation to coordinate space and back. The result is 

T{a x + a 2 - d/2)T{d/2 - a x )T{d/2 - a 2 ) ... 

G{a u a 2 ) = =7 — ■-=-, — ^=-, r . (3) 

1 (ai ji (a 2 jl [a — ai — a 2 ) 

All one-loop integrals with integer indices can therefore be expressed in terms of 
G = T(l + e)r 2 (l - e)/r(l - 2e), with rational coefficients. 




Figure 1: Massless two-point diagrams 

Two-loop integrals of the form of Fig. [|b are evaluated by repeated use of (^) . 
The integral of Fig. [l]c is non-trivial. We write it as 

d d k d d l 1 

= F(ai,a 2 ,a 3 , a 4 , a 5 ). (4) 

Integration by parts gives the recurrence relation 

(d - ai - a 3 - 2a 5 )F = [ail+(5~ - 2") + a 3 3+(5~ - 4")] F, (5) 

where, for example, 1 + increases a±, and 2~ decreases a 2 . Applying this relation 
sufficiently many times, we can kill one of the lines 2, 4, or 5, in Fig. Qc. The 
resulting integrals, in Figs. 0d,e, are trivial: 

F(ai, a 2 , a 3 , a 4 , 0) = G{a\, a 3 )G(a 2 , a 4 ), 

F(ai, a 2 , a 3 , 0, a 5 ) = G(ai + a 2 + a 5 - d/2, a 3 )G(a 2 , a 5 ). (6) 

All two-loop integrals can be thus expressed in terms two structures: G 2 , and 
Gi = r(l + 2e)r 3 (l - e)/rLl - 3e), with rational coefficients. This algorithm! 
is implemented in REDUCE0I, and as part of the 3-loop package MINCER0. 

3. HQET two point integrals 

We write the HQET one-loop integral of Fig. ^a, with arbitrary indices, as 
f d d k 1 f _u y 2 _ i(-2m) d - 2 ^ 

J (2Tr) d (-fc 2 )«i \k + Lu) (47r)<V 2 ( - a ^ a2 >- VI 
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It can be calculated, using the parameterization 

oo 

1 T(a + 0) f yP-Wy 



a a bP T(a)T(/3) J (a + by) a +P 1 
o 

where y has dimensions of energy. The result is 

r(2«i + o 2 - d)T(d/2 - ai) 



(8) 



I (ai, 02) = 



r(ai)r(a 2 ) 



(9) 



All one-loop integrals with integer indices can be expressed in terms of Iq = T(l 
2e)r(l — e) with rational coefficients. 






\a 2 






«3 


a 4 



Figure 2: Two-point diagrams in HQET or on-shell QCD 

The two-loop integrals of Figs. |]b,c are trivial. Multiplying the integrand of 
Fig. |d by 

1 

1 = - Oo + w) + (l + w) - (k + lo+ w) , (10) 

we can kill one of the heavy lines, reducing this integral to trivial ones, like those 
of Figs. |2|f,g. Only the integral of Fig. ||e is non-trivial. We write it as 



d d k d d l 



1 



(27r) d (27r) d (-fc 2 ) Q i(-/ 2 )^(-(fc_/)2)a 5 \k + uj J \l + u 



(47r) d 

The recurrence relation 



-1(01,02,03,04,05). 



(11) 



(d — ai — 03 — 04 — 2as + 1)1 

= [ail + (B~ - 2~) + (2(d - ai - a 2 - a 5 ) - o 3 - o 4 + 1)4~] J (12) 

allows us to kill one of the lines 2, 4, or 5, in Fig. ||e. The resulting integrals, in 
Figs. gf,g,h, are trivial: 

1(01,02,03, "4, 0) = 7(oi, o 3 )/(a 2 , 04), 

I(ai, 0, a 3 , 04, ag) = I(ai, a 3 + a 4 + 2o 5 - d)/(a5, a 4 ), (13) 
1(04,02,03, 0, a 5 ) = 7(ai + a 2 + o 5 - d/2, a 3 )G(a 2 , a 5 ). 



4 Multiloop Calculations in Heavy Quark Effective Theory 



All two-loop integrals can thus be expressed in terms of Iq and I\ — T(l + 4e)r 2 (l — 
e), with rational coefficients. 

This algorithm was proposed inO, where it was implemented^ in REDUCE. It 
was used for calculationof: the heavy-field anomalous dimensions (previously con- 
sidered in Wilson-lincEa and on-shellQ QCD calculations); the heavy-lighL cur- 
rent anomalous dimensionO (also obtained using Feynman parameterizationEJ and 
Gegenbauer polynomial expansion, with a more difficult momentum routingo) ; and 
the correlator of two heavy-light currentsEHIll. The coefficient of a s /n in this cor- 
relator is of order 10, making sum rules for heavy-light mesons unreliable, in our 
opinionEj. 

4. HQET three point integrals expanded in v' ~ v 

Here we consider three-point diagrams, containing HQET propagators with two 
velocities, v and v' , and depending on two residual energies, uj and uj' . Their 
expansions in v' — v to any finite order can be evaluated algebraically. The heavy- 
quark propagators are expanded in v' — v, and the resulting tensor integrals are 
reduced to scalar integrals. The one-loop integral of Fig. with arbitrary indices 
specifying the dependence on lo and uj' , is 

f d d k 1 / uj y 2 ( lo 1 \" 3 _ i{-2uj') d - 2ai T ( W\ 
J (2ir) d (-k 2 )^ \k^ToJ~) {k^+ly) (4^) d / 2 / (, a i' a 2,a 3 ;^ 7 J , 

J(oi, 02,03; x) = I(oti, «2 + ct3)x a2 2-Fi(2ai + a 2 + 03 - d, a 2 ;a 2 + 03; 1 — x). 

We write the two-loop integral of Fig. ||e as 
d d k d d l 1 



(27r) d {2n) d (-k 2 ) a ^(-l 2 ) a ^(-(k - l) 2 ) a -- \k + uj J \l + uj' 

^2uj') 2{d - ai - a2 - a ^ ( uj\ 
(4^? I(ai,a 2 ,a 3 ,ai,a 5 ;—). 



(14) 



In the Feynman gauge, all diagrams generate scalar integrals with 01.2, 5 < 1. More- 
over, the cases with 0:3.4 > 1 can be obtained from I(ai, 02, 1, 1, 05; x) by differen- 
tiation with respect to uj and uj' . Then, either some of the lines are already absent, 
or one of the lines 1, 2, 4, 5 can be killed by applying the recurrence relationQ 



05 - «i H (2«i + o 5 + o 3 — d) 



LO 
J 



o 5 ( 1 - ^ ) 5+(l" - 2") + ail + (2- - 5-) + — a 3 3+4- 



I. (15) 



This produces the trivial integrals of Figs. [|f,h, and the integral of Fig. [2^, given 
by 7(05,04)1(01, 03,04 + 205 -<L§r)- 

This algorithm was proposed inta; we are not aware of a package that implements 



it. The method was applied for checkingO, at small velocity transfer, the gejaeral 
result for the heavy-heavy current anomalous dimensionlla, and for obtainingu the 
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correlator of two heavy-light currents and one heavy— heavy current, at small v' — v 
(subsequently obtained, for all velocity transfersEJj, using methods based on differ- 
ential equations). Again there are large radiative corrections. However, these are 
largely cancelled in the ratio of the three-point and two-point sum rules, making 



,hrec 



predictions for the Isgur-Wise functioroEla relatively safe. 
5. On shell massive two point integrals 

Expansions of massive two-point diagrams to any finite order in p 2 — m 2 can be eval- 
uated algebraically. Such on-shell diagrams are necessary for QCD/HQET match- 
ing. We write the one-loop integral of Fig. |]a, with arbitrary indices, as 

d d k 1 im d-2( ai+a2 ) 

-M(a u a 2 ). (16) 



{2n) d (~k 2 ) a i(-(k + mv) 2 +m 2 ) Q 2 (47r) d / 2 

It can be evaluated by Feynman parameterization. Alternatively, we can use the 
inversion fc M — » k^/k 2 , after transforming to a dimensionless euclidean loop momen- 
tum. This transforms the on-shell massive denominator into an HQET denomina- 
tor, and we obtain 

.,, , , r(ai + a 2 - d/2)T(d - 2ai - a 2 ) , , 

M(ai,a 2 ) =I(d-ai -a 2 ,a 2 ) = — — — — r . 17 

1 (a 2 )l (a — ai - a 2 ) 

All one-loop integrals with integer indices can be expressed in terms of Mq = r(l+e) 
with rational coefficients. 

There are two distinct non-trivial topologies for two-loop integrals. We shall 
start with the simpler type M, of Fig. |^e, 

d d k d d l 1 



/ 



{2n) d (2Tr) d (-fc 2 )«i {-l 2 ) a i(-(k - l) 2 ) a s(-(k + mv) 2 + m 2 )^ 

1 m 2(d -I>) 

x ( ii i Y2~, 2^7 = 7T~Td — *f (ai, a 2 , «3, «4, as). (18) 

(— (t + mv) 1 + m z ) ai (47r) a 

Using inversion, we can relate it to an HQET two-loop integral 

M(ai,a 2 , as, ot±, as) = I(d — ai — «3 — «5, d — a 2 — on — as, as, on, as). (19) 

However, this HQET integral contains two non-integer indices, and is hence more 
difficult to calculate than before. The recurrence relation for M has a form identical 
to (||), allowing us to kill one of the lines 2, 4, or 5 of Fig. ^e. The integrals of 
Figs. ^f,h are trivial, as in (fl3|), but that of Fig. ||g is not. Special combinations of 
recurrence relations allow us to kill one more line. As a result, all two-loop integrals 
of the type M can be expressed in terms of M5 2 and Mi — T(l + e)r 2 (l — e)T(l + 
2e)r(l - 4e)/(r(l - 2e)r(l - 3e)) with rational coefficients. 

We write the most difficult integrals, of the type N of Fig. ||d, as 

d d k d d l 1 

(27r) d (2?r) d (-k 2 ) ai (-l 2 ) a ^ (—(k + mv) 2 + m 2 ) a 3 (-(/ + mv) 2 + m 2 )"" 

1 m 2 ( d -E«) 

/ / t i t i = 7a \d — A(ai,a2,a 3 ,a4,a 5 ). (20) 

(— [k + I + mv) 1 + m z ) a& (47r) a 
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The recurrence relation 

(d - 2a 2 -oi4- a 5 )N = [a 4 4+2- + a 5 5+(2" - 3~)]N (21) 

allows us to kill one of the lines 2 or 3 of Fig. ||d. This gives integrals shown in 
Fig. ||g, belonging to the class M, and integrals shown in Fig. ^j, which are still non- 
trivial. Special combinations of recurrence relations allow us to kill one more line. 
As a result, all two-loop integrals of the type N can be expressed in terms of Mq , 
Mi, and N(0, 0, 1, 1, 1) with rational coefficients. Instead of using A^fLO, 1, 1, 1) as 
a basis integral, it is more convenient to use the convergent integralE3a 

N(l, 1, 1, 1, 1) = tt 2 log2 - |C(3) + 0(e), (22) 

whose 0(e) and 0(e 2 ) terms have been found by algebraic methods!, but are not 
necessary for two-loop calculations. 

This algorithm is implemented in the REDUCE package RECURSORej and the 
FORM package SHELL2!. It was used for calculating heavy-quark on-shell massi 
and wave-functioni renormalizations, and the QCD/HQET matching of heavy- 
light currents!. Coefficients of (a s /7r) 2 in the matching coefficients are typically of 
order 10, leading to substantial two-loop corrections in ratios such as /g» //b- 
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